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The D* D* and B* B* systems are studied dynamically in the one boson exchange model, where 
n, n, a, p and ui exchanges are taken into account. Ten allowed states with low spin parity are 
considered. We suggest that the 1 , 2 ++ , ++ and _+ B*B* molecules should exist, and the 
D* D* bound states with the same quantum numbers very likely exist as well. However, the CP 
exotic + , 2 + ~) B* B* and D* D* states may not be bound by the one boson exchange potential. 
We find that the 1 = configuration is more deeply bound than the 1 = 1 configuration, hence 
Z^(4050) may not be a D*D* molecule. Although y(4008) is close to the D* D* threshold, the 
interpretation of Y (4008) as a D* D* molecule is not favored by its huge width. 1 D* D* and B* B* 
states can be produced copiously in e + e~ annihilation, detailed scanning of the e + e _ annihilation 
data near the D* D* and B* B* threshold is an important check to our predictions. 

PACS numbers: 12.39.Pn, 12.39.Jh, 12.40.Yx, 13.75.Lb 



I. INTRODUCTION 

In the past years, the observations of a number of charmonium-like " X, Y, Z" mesons at B factories have stimulated 
the interest in the spectroscopy of the charmonium states again. There is growing evidence that at least some of these 
new states are non-conventional cc states, such as deuteron like hadronic molecules, tetraquark states or hybrid have 
been suggested Among these new mesons, some are very close to the threshold of two charmed mesons, such as 
X (3872) H and Z+(4430) 0,11 . This distinctive character inspires the molecular interpretation for these mesons. In 
particular, some new enhancements near the D*D* threshold have been observed recently. 

The Belle collaboration reported a broad ir + it~ J j 'ip peak near 4008 MeV in addition to the well-known state 
y(4260) by studying the initial state radiation process e + e~ — > 7/s.R7r + 7r~ J/tp [f|, its mass and width are fitted to 
be M = (4008± 40t28 4 ) MeV and r = (226 ± 44 ± 87) MeV respectively. We notice its width is huge. This peak was 
suggested to be related to the D*D* threshold and could be a D*D* molecule in Ref. [||. The 7r + 7r~ J/V> spectrum 
was studied further by the Babar collaboration. However, there was no evidence for this broad enhancement, and the 
an upper limit B(tt + tt~ J/ip)T e + e - < 0.7eV at 90% C.L. was obtained 0. F(4008) is far from being established so 
far, more experimental efforts are obviously needed. 

Of special importance is the observation of the state carrying non-zero electric charge with hidden charm quarks, 
bmcc the observation of Z+(4430) 0Q, the Belle collaboration reported two resonance-like structures Z^~ (4050) and 
Z^~(4250) in the ir + Xci mass distribution in the exclusive process B° — ► K~n + Xci §• Their masses and widths are 
determined to be Mi = (4051 ± 14^?) MeV, Ti = 82±^±g MeV, M 2 = (4248t^t|f°) MeV and T 2 = (177±||±|i 6 ) 
MeV respectively Since tt + is an isovector with negative G-parity, and \ci is a isospin singlet with positive G-parity, 
the quantum numbers of both Z^~(4051) and Z \ (4250) are I G = 1~. In Ref.0, Z^ (4050) was suggested to be possibly 
a J p = + D*D* molecule due to its closeness to the D*D* threshold. However, the QCD sum rule results indicated 
that the D* + D*° state is probably a virtual state, which is not related with the .2^(4050) resonance-like structure 
[Tol ] . In addition, we demonstrated that Z^"(4250) as a D\D / ' DD\ or DqD* / D*Dq molecule is disfavored in Ref. [TT| | . 

Since the repulsive kinetic energy is greatly reduced by the heavy quark mass, the interaction between light quarks 
is strong enough so that the molecular states consisting of heavy flavor mesons very likely exist. In fact, the hadronic 
molecules consisting of two charm mesons were suggested long ago [T3|. De Rujula, Georgi and Glashow proposed 
that the molecular states involving hidden cc pair do exist, and have a rich spectrum [l3[. Possible new resonance 
near the D*D* threshold was suggested by Voloshin fLU fl5|. Different from other possible exotic structures, there 
are uncontroversial evidences for hadronic molecule such as the deuteron, which is unambiguously a proton-neutron 
bound state. The deuteron has been studied in great details over the years 0, From these studies, we learn that 
the pion exchange determines most of the binding energy and the long range part of the deuteron wavefunction, and 
the S-D wave mixing effect plays a critical role in providing the binding. Guided by the binding of deuteron, Tornqvist 
performed a systematic study of possible deuteronlike two mesons bound states with long distance one pion exchange 
[ill, Oil . At short distance, the interaction should be induced by the interactions among the quarks. However, a detailed 
and reliable modelling of the short range interaction is not a easy matter, and various phenomcnological models have 
been proposed [20L [2l1|. although one pion exchange is expected to be dominant for the hadronic molecule. Inspired 
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by the nucleon-nucleon interactions, we further extended the one pion exchange model to include the short distance 
contributions from the heavier bosons 77, <7, p and uj exchanges in Ref. [22| . We have also taken into account the 
contribution of the " 8 function" term, which leads to a 8 function term in the effective potential in configuration space 
when no regularization is used. This one boson exchange model gives a very good description of the weakly bound 
hadronic molecule. It has been successfully applied to exploring the possible heavy flavor pseudoscalar- vector molecular 
states [22|, the D*D* system and the molecular interpretation of y(4140) Motivated by the controversial states 
y(4008) and (4050), we shall examine for which quantum numbers the boson exchange potential is attractive 
and strong enough so that bound states are expected, ten allowed D* D* states with low spin parity are considered. 
Moreover, the B*B* system would be discussed as well. 

The paper is organized as follows. In section II, the formalism of one boson exchange model is summarized. In 
section III, we apply the one boson exchange model to the D*D* system, the quantum numbers of the D*D* bound 
states which might exist, are suggested. The B*B* system is discussed along the same line in section IV. Finally we 
present our conclusions and some discussions in section V. 



II. THE FORMALISM OF THE ONE BOSON EXCHANGE MODEL 



In the one boson exchange model, the effective potential between two hadrons is obtained by summing the in- 
teractions between light quarks or antiquarks via one boson exchange. To leading order in the boson fields and 
their derivative, the effective interactions between the constituent quark and the exchanged boson are as follows 

[mums! 



Pseudoscalar : C p = — g pqq ip(x)ij5i^(x)ip(x) 

Scalar : C s = — g sqq ^{x)^(x)4>(x) 

Vector : £ v = -g vqq ^(xh^(x)v»(x) - ^^(x)a^(x)dV (x) 



Here m q is the constituent quark mass, ip{ x ) is the constituent quark Dirac spinor field, (p(x), 4>{x) and v fJ, (x) are 
isospin-singlet pseudoscalar, scalar and vector boson fields respectively. In this work we take m q = m u = ~ 313 
MeV, since we concentrate on the constituent up and down quarks here. If the isovector bosons are involved, the 
couplings enter in the form t • if, t • <f> and r • v 11 respectively, where r is the well-known Pauli matrices. Following 
standard procedure, we can straightforwardly obtain the one boson exchange potential between two quarks. 

1. Pseudoscalar boson exchange 

2 3 

Vp ( T ) = 9 ™ i [ _ Hl (A, m p ,fx p , r) <r, ■ a, + H S (A, m p ,fi p , r)S« (f )] (2) 

where SV,(f) = 3(<r.; • r)(<Xj • f) — cr; • er 3 is the tensor operator. We have defined fj, p = m p — (myi — mv2) 2 to 
approximately account for the recoil effect due to the small mass splitting within the heavy flavor vector meson 
isospin multiplet (l8l . Il9l . [27j . where m p is the exchanged pseudoscalar mass, myi and my2 are the masses of 
the heavy flavor vector mesons involved. For the B*B* system, the mass splitting can be negligible so that 
m„ is satisfied. 



2. Scalar boson exchange 



9 sqq 



p 2 /i 2 

H Q (k,m s ,p, s ,r) + -^-^Hi(k,m s ,ii s ,r) + -^-^H 2 (k, m s , /i s , r)L • Sy 



where Sy = + erj), p? s = m 2 — (myi — mv2) 2 with m s being the exchanged scalar meson mass, and 

L = r x p is the angular momentum operator. 

3. Vector boson exchange 

{2 
9l qq Ho{K,m v ,p v ,r) - {gl qq + 4g vqq f vqq )-^Hi(A,m v , fj, v ,r) 

-(9vqq + f v <i<)) 2 ~[^2 #3(A,m„, fj, v , r) Sij (r) + 2H 1 (A,m v ,fi v ,r)(cr i ■ a 

(4) 
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where (j% — ra 2 v ~~(mvi—niv2) 2 approximately reflects the recoil effect with m v being the exchanged vector meson 
mass. For 1=1 isovector boson exchange, the above three potentials in Eq. (J2J)-Eq. should be multiplied by 
the operator Tj • Tj in the isospin space. 

The dimensionless functions H (A, to, fi, r), -ffi(A, to, p, r), ^(A, m, p,r) and H 3 (A, to, p, r) introduced in Eq.([2])- 
Eq.(|D) are defined as follows 

Ho(A,m, M ,r) = - (e~^ - e~ Xr ) - ^f£^ e~ Xr 
pr v ' 2pX 

Hi (A, to, h, r) = -l-{e-^-e-^)+ X{K2 - m2) e~ x * 
pr 2p A 

2 ™2 .-Xr 



1 \ 1 ,„ 1 J 1 _ Xr A — TO e 



F 2 (A,TO )M ,r) = (l + _)_ e -^-(i + _)__ e 



pr p 2 r 2 /i /z 2 r 2 2/i 2 /ir 

H 3 (A,TO )M ,r) = + 3 } ^ e _ Mr _ (l+ 3 + 3 ^e^_A 2 _^ 

yur p z r z pr Xr X z r z p z /ir 

with AT 2 = A 2 + p 2 — to 2 . In deriving the above effective potentials, we have introduced form factor at each interaction 
vertex to regularize the effective potential at short distance, and the form factor in momentum space is taken as 

m = (6) 

where A is the so-called regularization parameter, to and q are the mass and the four momentum of the exchanged 
boson respectively. This form factor suppresses the contribution of high momentum, i.e. small distance. The presence 
of such a form factor is dictated by the extended structure of the hadrons. The parameter A, which governs the 
range of suppression, can be directly related to the hadron size which is approximately proportional to 1/A. However, 
since the question of hadron size is still very much open, the value of A is poorly known phenomcnologically, and it 
is dependent on the models and applications. In the nucleon-nucleon interactions, the A in the range of 0.8-1.5 GeV 
has been used to fit the data. For the present application to the heavy flavor vector mesons system, which have a 
smaller size than the nuclcon, we would expect a larger regularization parameter A. We have demonstrated that the 
binding energy and static properties of the deuteron are produced very well in the one boson exchange model, if A is 
chosen to be about 808 MeV [22j . In the case that all coupling constants except g^NN are reduced by half, A should 
be approximately 970 MeV. The extended structure of hadrons also has the following obvious consequence: because 
the mass of the exchanged meson determines the range of the corresponding contribution to the D* D* interactions, 
one should restrict oneself to meson exchange with the exchanged meson mass below a certain value, typically on the 
order of the regularization parameter A. Since n, rj, a, p and u> exchanges are considered in the present work, the 
value of A should be larger than the u meson mass. 

In the present one boson exchange model, the input parameters include the masses of the exchanged bosons and 
heavy flavor vector mesons, and the effective coupling constants between the constituent quarks and the exchanged 
bosons. The meson masses are taken from the compilation of the Particle Data Group |28j: m 7T ± = 139.57 MeV, 
m w o = 134.98 MeV, m v = 547.85 MeV, m a = 600 MeV, m p = 775.49 MeV, to w = 782.65 MeV, m D ,o = 2006.97 MeV, 
m D ,± = 2010.27 MeV and tob* = 5325.1 MeV. The constituent quark-meson coupling constants can be estimated 
from the phenomcnologically known ttNN, rjNN, <tNN, pNN and ojNN coupling constants via the well-known 
Goldberger-Treiman relation [22|, [29[ . 

3 TTlq rriq 

g-xqq = 7 SttAW, g-nqq = g v NN 

5 TOJV TOjV 

3 Tfb 3 TYl 

9pqq = 9pNN, fpqq = - —^fpNN - (1 - ~— —)g P NN 

O TOJV O TOjv 



1 „ m„ „ ,1 to 

3 
1 



„ gcuNN, fujqq = -JujNN ~ (~ ~)guiNN 

6 tun o m,M 



daqq ~ g g<jNN (7) 

Among the above effective boson- nucleon coupling constants, only g^NN has been determined accurately from the 
pion-nucleon and nucleon-nucleon scatterings. As in our previous work [22I l23l|, the effective coupling constants are 
taken from the famous Bonn model pr| . The uncertainty of the coupling constants will be taken into account later, all 
the coupling constants except g^NN are reduced by half for demonstration, and the corresponding numerical results 
arc presented as well. 
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D* D* 



Tr,ri,a, p,tv 



D* D* 

FIG. 1: D* D* interactions in the one boson exchange model at quark level, where the thick line represents heavy quark or 
antiquark, and the thin line denotes light quark or antiquark. 

For the system consisting of two vector mesons, the spatial parity is determined by P = ( — 1) L and the C parity is 
C = ( — 1) L+S , where L is the relative angular momentum between two vector mesons, and S is the total spin of the 
system. We cutoff the total angular momentum of the system to J = 2, the allowed states with low spin parity are 
listed in Table [U In the following, we shall study for which quantum numbers the one boson exchange potential is so 
attractive that bound states may exist. 





Channels 


0++ 


i 5o, a Do 


1+- 


3 S 1 , 3 D 1 


o-+ 


3 Po 


1++ 


5 £i 


1-+ 


3 Pi 


2+- 


3 D 2 


1"- 


1 P L , 5 Pi, B Fi 


2++ 


1 D 2 , 5 S2, ''D2, °G2 


2"+ 


3 P 2 , 3 F 2 


2" 


5 P 2 , 5 F 2 



TABLE I: The allowed states of the system consisting of two vector mesons, where only states whose total angular momentum 
is smaller than 3 are listed. 



III. THE POSSIBLE MOLECULAR STATES OF THE D*D* SYSTEM 

There is a sign difference (— I) G between the quark-quark interactions and quark-antiquark interactions, the mag- 
nitudes are the same, where G is the G-parity of the exchanged meson. Consequently both 7r and oj exchanges give 
opposite sign between quark-quark interactions and quark-antiquark interactions. In the present case, the effective 
potential is induced by one boson exchange between a pair of light quark and antiquark, and the diagram contributing 
to D*D* interactions is shown in Fig. [TJ The effective potential is explicitly expressed as 

V(t) = -V^T) + V v {v) + V a {r) + V p {y)-V^v) 

= Vc(r) + Vs{r){(Ti ■ o-j) + Vi(n p ,r)(Ti ■ tj) + V T (r)Sij(r) + Fs/(m w , fj, n , n p ,r)(tTi ■ <Tj)(Tj • t 3 ) 

+V T i(m 7r , fj, n , n p ,r)Sij(r)(Ti ■ tj) + V LS {r)(L ■ Sy) + V LS i(fJ- P ,r)(L ■ S^)^ • t 3 ) (8) 

where Vm i r ) {■M = it, r), cr, p and uj) denotes the effective potential induced by the meson M. exchange between two 
quarks. The subscripts i and j are the indexes of light quark and antiquark. The spin operator (isospin operator) rr^ 
or <Tj (xj or Tj) only acts on the light quark and antiquark. The parameters m w , p^ and fi p could take two different 
sets of values due to the small mass splitting within the D* and w isospin multiplets. For ]J* D* — > ]J* D* and 
£)*+D*~ — > D* + D*~, we should choose m T = m n o, fa = ra„n and p p = m p . Whereas for the D*°D*° — > D* + D*~ 
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and D*+D* 



0*°D* processes, we should take = m^±, fi v = [r 



+ — mu«o 



,211/2 



and 



211/2 = 



A*p = [ m l - (mD*+ - m D *o) 



Vc(r) = 



fjL Pl . The eight potential functions Vc(r), Vs(r) etc are given by 



Hq(A, m a ,m a ,r) 



8m 2 . 



H 1 (A 1 m a ,m IJ ,r) 



in 



■m w H {A,m u ,m w ,r) 



9ujqq ^dujqqfujqq 771^ 

47T 8m 2 , 



.ffi(A, m u , m u ,r) 



Vs(r) 



Sijqq m rj 

^bTl2m 2 



Hi(A,m v ,m v ,r) 



{9^qq "F fuqq) m u 



in 



6m 2 



Vj^r) = 9 -^-^ p H (A,m p ,(x p ,r) - 



+ 



pqqfpqq /^p , . . 



V r (r) 



4^ 12TO 2^3(A,m 1) ,m 1) ,rj + 



4tt 

"I" fujqq) 



in 



I2m 2 q 



H 3 (A, m^,m a 



g 2 3 

Vsi{m„, 11^,11 p ,r) = n Hi (A, , , r) 

47r 12m~ 



(fiWg + fpqq) 2 M, 



4tt 



6m 2 



—Hx(A,m p ,n P ,r) 



9nqq ^4 u / A „, ,, (ffpgg + f pqqY 



MP 



4tt 



12m 2 



H 3 (A,m p ,ti p ,r) 



V LS (r) 



?ggg 



in 2m 2 



Z-H 2 (A, m a ,m 



3.9, 



2 

"gg 



^Quiqqfu. 



V L si(p P ,r) = - 



^9pqq 



^9pqqfp 



in 



2m 2 , 



in 



H 2 (A,m p ,fj, p ,r) 



2m 2 



H 2 {A,m UJ ,m u>1 r) 



(9) 



Since the threshold of D* + D*~ is about 6.6 MeV higher than the D*°D*° threshold, the isospin symmetr y is expected 
to be violated drastically for the D*D* molecular states whose binding energy is of order a few MeV [22|]. Isospin 
violation mainly comes from three aspects: the first is the different kinetic energies for £)*°£)* and D* + D*~, the 
second is the different effective potentials from tt° exchange and tt^ exchange, and the third is because of the different 
thresholds of ]j* f)* and D* + D*~ . In the following, we will perform the same analysis as that for the deuteron and 
the possible heavy flavor molecules in Ref. [13, HJ. One can then determine for which quantum numbers the one 
boson exchange potential is attractive and strong enough so that the D*D* bound states are expected. Firstly we 
consider the J PC = 0++ D*D* states as an demonstration, the system can be in S wave or D wave similar to the 
deuteron. Taking into account the isospin violation effect, one has four coupled channels. For convenience, we choose 
the basis to be |1) = I 1 S a (D*° 'D* )} , |2) = \ 5 D (D*°D* )), |3) = \ 1 S (D*+D*-)) and |4) = \ 5 D a (D*+D*-)), then 
the wavefunction of the system is written as 



\0 ++ (D*D*)) 



\ L S {D* U D*")) + 



u 2 {r) 



\ L S Q (D*+D*-)) + 



^1\ 5 D (D*+D*-)) (10) 



where u±(r), u 2 (r), 1*3(7") and U4 (r) are the spatial wavefunctions. The matrix elements of the light quark relevant 
operators <x, -<rj, Sij (r) and L-S^ etc in Eq.® can be calculated straightforwardly with the help of angular momentum 
algebra, and the results are given analytically in the Appendix of Ref. [S^]. Consequently the one boson exchange 
potential for the ++ D*D* state can be written in the matrix form as 




V s {r) 



( -2 
10 
0-20 

V 1 



Vi(/J, P ,r) 



1 





-2 








-1 





-2 


-2 





-1 








-2 





-1 



-V T i{m^,^^ p ,r) 






-V2 
-2 / 

/ V2 2^2\ 
V2 2 2^2 4 

2V2 V2 
V 2V2 4 V2 2 ) 



V LS (r) 







/ 2 4 
0-10-2 
4 2 

\0 -2 -1 


0-300 



-3 



+ V LS i{n P ,r) 



/o 
3 6 


\ 6 3 



(11) 
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For the up- left and down- right 2x2 matrix elements, we should choose m n — m w o, fi n — m v o and /i p = m p . While for 
the off-diagonal 2x2 matrix elements, we should take m n = m 7T ± , fj,^ = fi 7T1 and fi p = /i pi . Taking into account the 
D wave centrifugal barrier and solving the coupled channel Schrodingcr equation numerically, the numerical results 
are listed in the Table |TTJ It is obvious that the binding energy and the static properties are rather sensitive to the 
regularization parameter A and the effective coupling constants, this is common to the one boson exchange model 
[22, HH, H3] ■ We also find that the binding energy increases with A, this is because increasing A increases the strength 
of the potential at short distance. For A = 930 MeV, a bound state with mass about 4013.80 MeV appears. We 
can see that the isospin symmetry is strongly broken especially for the states near the D*D* threshold. Fig. [5] 
displays the wavefunction of the bound state with mass 4004.40 MeV and A = 950 MeV. One notice that the D*°D*° 
component dominates over the D* + D*~ component for both the S wave and D wave configurations, as could be 
expected. Because the wavefunctions U\{r) and u s (r) have the same sign, the same is true for U2{r) and u^(r). 
Therefore the 1 = component is predominant for this state, it would be an isospin singlet in the isospin symmetry 
limit. The dominance of the 1 = configuration is observed for all the states listed in Table [TT1 From the numerical 
results, we can see that the D wave probability increases with the regularization parameter A, the importance of 
the tensor force is obvious. The uncertainties induced by the effective coupling constants are considered as well. All 
the coupling constants except g^NN are reduced by half, and the corresponding numerical results are presented in 
Table |TTJ The same pattern of the static properties dependence on A is found. Bound state solution appears if the 
regularization parameter A is about 1100 MeV, the value of A is still in the reasonable range. Since the molecular 
state is widely extended, the decay into light mesons via annihilation is generally suppressed by the form factor. The 
leading source of decay is dissociation, to a good approximation the dissociation will proceed via the almost free space 
decay of the constituent mesons. Consequently the ++ D*D* molecule mainly decays into DD-f-f, and DDjit, and 
the mode DDtttt is strongly suppressed or even forbidden by the phase space. 




r(fm) 



FIG. 2: The four components spatial wavefunctions of the ++ D* D* state with A = 940 MeV. 

The one boson exchange potentials for the nine remaining states are listed in Appendix A. Following exactly the 
same method, the binding energy and the static properties can be predicted, and partial results are shown in Appendix 
B for illustration. The binding energy and the static properties are found to be rather sensitive the regularization 
parameter A. The spatial wavefunctions for the ]J*°£)* and D* + D*~ components have the same sign, consequently 
the 1 = component is dominant for all these states, and they are isospin singlets in the isospin symmetry limit. The 
same conclusion has been reached in the one pion exchange model where the strength of the effective potential 
for the 1 = state is about one third of that for / = 1. We shall discuss these states one by one in the following. 

For the axial vector 1+" state, there are four channels 3 S 1 (D*°D*°), 3 D 1 {D*°D*°), 3 S 1 {D* + D*~) and 
3 Di(D* + D*~). The energy of the system is substantially lowered due to the S-D wave mixing effect. The cou- 
pling between the S wave and D wave has the same strength as the ++ case, which is clearly seen from Eq. lfTTj) and 
Eq. (|Al[) . It is obvious that the predictions for the binding energy and the static properties have similar pattern with 
the ones for the ++ state, and the binding of the 1 H state is less stronger than the ++ one for the same A value. 
For A as large as 980 MeV, we can find a bound state with mass about 4012.43 MeV. If all coupling constants except 
9ttNN are reduced by half, bound state begins to appear for A ~ 1200 MeV. We note that the unnatural spin parity 
forbids its decay into DD, while the decay mode DD* / D* D is allowed. 

The pseudoscalar h D* D* state involve two channels 3 Pq(D*°D*°) and 3 P (D* + D*~). If the small isospin 
violation effect is neglected, the two coupled channel problem is reduced to a familiar one channel problem. Although 
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there is repulsive P wave centrifugal barrier, the one boson exchange potential is so strong that the P wave centrifugal 
barrier can be partly compensated, therefore bound state solutions can be found for reasonable values of A, as can be 
seen from Table ITTll For A = 950 — 1030 McV, we find that the binding energy with respect to the D*°D*° threshold 
is in the range of 2.6 to 137.2 MeV. The binding energy is more sensitive to A than the ++ and 1^ four coupled 
channels cases. The probabilities for the 3 Po(D*°D*°) and 3 Pq(D* + D*~) components are close to each other, as 
could be expected. The small difference is induced by the mass splitting within the D* and tt isospin multiplet. 

The results for the 1++ state are similar to the h case. Because the D wave centrifugal barrier is higher than 
the P wave centrifugal barrier, the total effective potential for the 1 ++ state is less attractive than the h state. If 
the coupling constants except g^NN are reduced by half, bound state solutions can be found only for A larger than 
1290 MeV. Since the numerical results indicate that larger value of A is required to bind the 1 ++ D*D* state, the 
1++ D*D* state is harder to be bound than the previous states considered. 

For the CP exotic 1 h and 2 H states, the system is in P wave and D wave respectively. Considering isospin 
violation effect, two channels are involved for both states. As has been shown in Eq. (|A4[) and Eq. (|A5[) . the one boson 
exchange potentials for these two states are exactly the same, and they are less attractive than the potentials for the 
h and 1 ++ states. For the 1 h state, bound state solution appears only for A as large as 1300 MeV, and we can 
find 2 H bound state only if the regularization parameter A is larger than 1640 MeV. If we reduce all the coupling 
constants except g^NN by half, A larger than 3010 MeV and 5110 MeV respectively for the 1 ^ and 2^ states is 
required to find bound state solutions. Because the value of A is so large that it is far beyond the range of 0.8 to 1.5 
GeV favored by the nucleon-nucleon interactions, we tend to conclude that the CP exotic 1 and 2 H D*D* states 
can not be bound by the one boson exchange potential. This conclusion is consistent with the fact that no such CP 
exotic states have been observed so far. 

We then come to the very interesting 1" D*D* state, there are six configurations 1 Pi(D*°D* a ), 5 Pi(D*°D*°), 
5 Fi(D*°D*°), 1 P 1 (D*+D* -), 5 P 1 (D*+D*-) and 5 F 1 ( D*+D*~). Due to the substantially strong attraction of the 
effective potential, bound state solutions can be found for reasonable value of A in spite of the P wave centrifugal 
barrier. The wavefunction of the bound state with mass about 4006.82 MeV and A = 920 MeV is displayed in Fig. 
[3J From the numerical results in Table HVl we see that the isospin symmetry is violated, especially for the states near 
the threshold. The 5 Pi(D*°D*°) and 5 Pi (D* + D*~) are the dominant components, the 5 Fi components are strongly 
suppressed by the large F wave centrifugal barrier. For the same value of the regularization parameter A, we notice 
that the binding energy of 1~~ state is the largest among the ten allowed D*D* states. Hence we suggest that the 1 _ ~ 
D*D* molecular state should exist, this conclusion is consistent with the results obtained from the general quantum 
mechanical properties of unitarity and analyticity [3, EE]- The 1 D*D* state is remarkable, it can be directly 
produced via the e + e~ annihilation or with the help of the initial state radiation (ISR) technique at B factory. The 
existence of such a state can be either confirmed or rejected if more detailed e + e~ annihilation data near the D*D* 
threshold become available. We strongly urge the Babar and Belle collaboration to search for this state. The 1~~ 
D*D* molecule mainly decays into DD-fj and DDjir via the dissociation of D* and D* , the decays into DD and 
DD* I D* D are allowed as well. The width of the 1 D*D* molecule should be of the same order as the D* width. 
Therefore it would be be a narrow state, and its width is expected to be of the order about 10 MeV. For the 1 _ ~ 
state y(4008) reported by the Belle collaboration, although it is close to the D*D* threshold, its width is huge, which 
is r = (226 ± 44 ± 87) McV. Consequently it seems unreasonable to interpret V(4008) as the 1 D*D* molecule. 
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FIG. 3: The spatial wavefunctions of the 1 D*D* state with A = 920 MeV. 
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For the 2++ state, eight channels 1 D 2 {D*°D*°), 5 S 2 (D*°D*°), 5 D 2 {D*°D*°), 5 G 2 (D*°D*°), 1 D 2 (D*+D*-), 
5 S 2 (D* + D*-), 5 D 2 ( D*+D*-) and 5 G 2 {D*+D*-) are involved. The effective potential in matrix form is given 
in Eq. (|A7j) . which is more complex than the previous cases considered. It is obvious that both the tensor interaction 
and the spin-orbit interaction vanish in the 5 S 2 configuration. However, bound state solution can be found for rea- 
sonable value of A (A should be larger than 860 MeV and 970 MeV respectively for the two sets of coupling constant 
values). The reason is that the mixing of 5 S 2 with 1 D 2 , 5 D 2 and 5 G 2 under the tensor force increases the binding 
of the system considerably through higher order iterative processes. We show the wavefunction of the 2 ++ molecular 
state with mass 4010.82 MeV and A = 860 MeV in Fig. 2] Isospin violation is obvious, and the I = component is 
dominant. From the numerical results presented in Table fVl we can see that 5 S 2 (D*° D*°) and 5 S 2 {D* + D*~) are the 
dominant components, the 5 D 2 probability is larger than the 1 D 2 probability, and the b G 2 component is suppressed. 
It is remarkable that the dependence of the 2 ++ state mass on A is the least sensitive among the ten allowed states. 
Comparing the numerical results presented in Table |TT] and Table [V] it can be seen that the binding of this state is 
comparable to the ++ state, accordingly the 2 ++ D*D* molecule very likely exist. 
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FIG. 4: The spatial wavefunctions of the 2 ++ D* D* state with mass 4010.82 MeV and A = 860 MeV. 

Finally two states with J PC = 2 h and 2~~ respectively remain. For both states, one has four coupled channels, 
P wave and F wave configurations are involved. Because of the P wave and F wave centrifugal barrier, bound state 
begins to appear for A as large as 1150 MeV and 1090 MeV respectively. If all the coupling constants except g W NN 
are reduced by half, the value of A should be larger than 1630 MeV and 1480 MeV respectively in order to find bound 
state solutions. We tend to conclude that the one boson exchange potential may not support the 2 h and 2~~ D*D* 
states. 

In short summary, ten allowed D*D* states with low spin parity has been studied. We find that the isospin 
symmetry is violated, especially for the states near the D*D* threshold. The 7 = configuration is dominant for all 
the lowest bound states, hence they would be isospin singlets in the isospin symmetry limit. Since (4050) is an 
isospin vector, 2^" ( 4050) as a D*D* molecule is not favored. This conclusion is consistent with the prediction from the 
QCD sum rule [HI. We suggest that the 1 — , 2++, 0++ and _+ D*D* molecules should very likely exist, whereas 
the CP exotic 1 h and 2 H D*D* states, 2 ^ and 2~~ D*D* states should not be bound by the one boson exchange 
potential. Although V(4008) is close to the D*D* threshold, its width is so large that it is not reasonable to identify 
V(4008) as a D*D* molecule. The 1 D*D* molecular state can be produced copiously in e + e _ annihilation or via 
the initial state radiation at B factory, detailed scan of the e + e~ annihilation data near the D*D* threshold is crucial 
to confirming this prediction. 

IV. THE MOLECULAR STATES OF THE B* B* SYSTEM 

The mass difference between B* + and B*° is so small that it can be negligible [28|. If we take into account the mass 
splitting within the exchanged pion isospin multiplet, we should solve similar coupled channel problems as the D* D* 
case. The repulsive kinetic energy is greatly reduced due to the larger mass of B* meson, therefore the B*B* system 
should be more deeply bound than the D*D* system. Numerically solving the corresponding Schrodingcr equation, 
we notice that the molecular states bound for reasonable A value have definite isospin, and the 1 = configuration 
is obviously much easier to be bound than the 1=1 configuration. We also find that the binding energy dependence 
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on A becomes less sensitive, if we reduce all the coupling constants except g^NN by half. In the case that the mass 
difference between tt ± and it is neglected, the B*B* state is of definite isospin, and the dimensions of the coupled 
channel equations would be reduced by half. We have seriously calculated the binding energy and the static properties 
for both the 1 = and 1 = 1 states. However, the numerical results are too lengthy to be listed in the manuscript. 
We find that introducing the pion mass splitting will modify the binding energy by at most 0.5 MeV. 

For the ++ and l" 1 B*B* system, both the 1 = and 1=1 bound states can be found for the same value of the 
regularization parameter A, whereas the two states behave in different way. The 1 = state is generally more deeply 
bound than the 1=1 state. For the isospin singlet, the D wave components increase drastically with A. Whereas 
for the isospin vector, the S wave components arc dominant, and they increase slightly with A. Because B* mainly 
decays into [lH, B.B77 is the leading decay mode of the B*B* molecule via the B* and B* dissociation. For the 
eight remaining states, we also find that the 1 = configuration is more tightly bound than the 7=1 configuration, 
hence it seems to be a universal result that the hadronic molecule prefers to a isospin singlet. 

Since the mass splitting within the pion isospin multiplct only introduces minor modifications. The numerical 
results can be understood easily in the exact isospin symmetry limit. If the pion mass splitting is neglected, the 
1 ++ , 1 h and 2 H B*B* states only involve one channel. Both the one boson exchange potential and the total 
potential including the centrifugal barrier for the four states are displayed in Fig. [5J From Fig. [3^, we can see that 
the P wave centrifugal barrier is partly compensated by the one boson exchange potential, and there remains a weak 
attractive interaction in the intermediate range. Therefore bound state can be found for reasonable values of 
A. For A = 808 — 950 MeV, we find the binding energy is in the range of 2.93 to 134.71 MeV. Because the D wave 
centrifugal barrier is higher than the P wave centrifugal barrier, the total potential of the 1 ++ state is less attractive 
than the ^ one for the same values of parameters, this point can be seen clearly by comparing Fig. 03 with Fig. 
[SJi. Accordingly our numerical results really indicate that the 1 ++ B*B* state is harder to be bound than the h 
state. For A = 900 MeV, 1 h B*B* bound state can not be found, it is because the potential including the P wave 
centrifugal barrier is repulsive in this case. When A is increased to about 1120 MeV, the total potential shown in 
Fig. E]i becomes attractive in the intermediate region, the 1 h B*B* system is marginally bound. The one boson 
exchange potentials for the CP exotic 1 h and 2 H states are exactly the same. However, the 2 H B*B* system is 
more loosely bound than the 1 h system because of the D wave centrifugal barrier. Fig. [5^ and Fig. clearly show 
that the total potential of the 2+" state is still repulsive for both A = 900 MeV and 1120 MeV. 2+" B*B* bound 
state begins to appear for A as large as 1260 MeV. 




02468 10 02468 10 2468 10 

r(GeV~') rfGeW 1 ) r(GeV~') 

(d) (e) (f) 

FIG. 5: The potentials for the single channel h , 1 ++ , 1 _+ and 2" 1 B*B* states with 1 = 0. (a) and (b) show the potentials 
for the h and 1 ++ states with A = 900 MeV respectively, (c) and (d) are the potentials of the 1 h states with A = 900 
MeV and 1120 MeV respectively, (e) and (f) are for the 2 + ~ states with A = 900 MeV and 1120 MeV respectively. The solid 
line represents the potential from one boson exchange, and the dashed line denotes the total potential including the centrifugal 
barrier. 



10 



The dominance of the 1 = configuration over 1=1 can be clearly understood as well for the h , 1 ++ , 1 h and 
2^ states in the isospin symmetry limit. In this case, the one boson exchange potential is the summation of the isospin 
independent part and Ci multiplying the isospin relevant part, where the parameter Cj is equal to -3 and 1 respectively 
for I = and 1. Concretely for the v state, the isospin independent potential is Vc{r) — Vs{r) — 2Vr{r) — V~Ls( r ), 
and the isospin relevant part is Vi(fi p , r) — Vgj(m w , fj^, fj, p ,r) — 2Vri(Tn x , fj, T , {i p , r) — Vlsi(Hp, r )- We plot both the 
isospin irrelevant and relevant potentials for these states in Fig. it is obvious that the isospin irrelevant potential 
is usually attractive. For the h and 1++ states with A = 900 MeV, the isospin relevant potential is positive, hence 
it is easily understood why the 1 = configuration is more attractive than 1=1. For the 1 h state with A = 900 
MeV, the isospin relevant part is negative, thus the potential for the 1=1 state is deeper than the one for 1 = 0. 
However, the corresponding potential still can not support a 1 h B*B* isovector state, such state can be bound only 
when A is increased to about 1200 MeV. For A as large as 1120 MeV, the isospin relevant potential becomes positive, 
accordingly the 1 = configuration is more tightly bound than 7=1. In this case, a marginally bound 1 h B*B* 
isospin singlet appears with mass about 10635.39 MeV. If we reduce all the coupling constants except g^NN by half, 
the 1 - + and 2+~ B*B* states can be bound only if A is larger than 1700 MeV and 2790 MeV respectively. Therefore 
the one boson exchange potential may not support the CP exotic 1 h and 2 H B*B* molecules. 




(a) (b) (c) (d) 

FIG. 6: The isospin relevant and irrelevant one boson exchange potentials for the h , 1 ++ and 1 h B*B* states, (a) and 
(b) respectively show the potentials for the ^ and 1 ++ states with A = 900 MeV. (c) and (d) are the potentials for the 1 _ + 
states with A = 900 MeV and 1120 MeV respectively. The one boson exchange potential for the 2 + ~ state is exactly the same 
as the 1 h case. The solid line represents the isospin irrelevant potential, and the dashed line denotes the isospin relevant part. 



Following the same methhod as the D*D* case and carefully examining the numerical results, we suggest that the 
1~~, 2 ++ , ++ and h B*B* molecules should exist, and the l" 1 B*B* bound state also very likely exists. Similar 
to the 1 D*D* state, the 1" B*B* molecule can be produced largely in e + e~ annihilation at Babar or Belle. 
Detailed scanning of the e + e~ annihilation data at the B*B* threshold is expected. 



V. CONCLUSIONS AND DISCUSSIONS 



Motivated by the charmonium-like state V(4008) and Z+(4050), the possible D*D* molecular states have been 
studied dynamically in the one boson exchange model, where x, ~q, a, p and uj exchanges are taken into account. Ten 
allowed states with low spin parity have been considered. We find that the binding energy and static properties are 
sensitive to the regularization parameter A and the effective coupling constants. The binding energy increases with 
A, whereas the root of mean square radius decreases with A, this is because increasing A increases the strength of 
the potential at short distance. Larger coupling constants are favorable to the formation of molecular states. If all 
the coupling constants except g n NN are reduced by half, larger value of A is required to find bound state solutions. 
Isospin violation is expected, especially for the states close to the threshold, and the 1 = component is dominant, 
The predominance of the 1 = configuration over 1 = 1 can be clearly understood in the exact isospin symmetry 
limit, and the same conclusion is reached in the one pion exchange model. Hence the interpretation of Z^~(4050) as 
a D*D* molecule is not favored. 

Since the regularization parameter A is poorly known so far, we are not be able to precisely predict the binding 
energies for the possible molecular states bound by one boson exchange potential. Certainly, if the potential is strong 
enough one can be quite confident that such bound state must exist, but their exact binding energy always depends 
on the details of the regularization procedure. However, we can reliably predict which ones of the ten allowed states 
are much easier to be bound, and the prediction is rather stable even if the uncertainty of the coupling constants is 
considered, as is obvious from the numerical results in Appendix [Bj Further research on X(3872) would put severe 
constraint on the parameters on the one boson exchange model, especially on the regularization parameter A, so that 
the predictions presented in the work could become more precise. 

Our detailed numerical results indicate that the l - ~, 2 ++ , ++ and 0~ + D*D* bound states should very likely 
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exist, whereas the CP exotic 1 h and 2 H D*D* states, 2 h and 2 __ D*D* states may not be bound by the one 
boson exchange potential. The 1 state can be directly produced in e + e~ annihilation, detailed e + e~ annihilation 
data near the D*D* threshold are important to confirm or refute the existence of such state. The D*D* molecule 
mainly decays into DD77 and DD-fir via the dissociation of D* and D*, the DDtttt mode is highly suppressed or 
forbidden by the phase space, and its width should be about tens of MeV. Although Y"(4008) is close to the D*D* 
threshold, its width is huge so that it is unreasonable to identify Y(4008) as a D*D* molecule. The existence of 
y(4008) and (4050) has not been confirmed so far, more experimental efforts are urgently needed. 

The possible B*B* molecular states have been discussed along the same line. The B*B* system is more deeply 
bound than the D*D* system. We find that the molecular states of the B*B* system have definite isospin, and the 
isospin singlet is much easier to be bound than the isospin vector for the states considered, the reason is analyzed 
in the isospin symmetry limit. If the small isospin mass splitting between the neutral and charged pion mesons is 
neglected, the dimension of the coupled channel problem is reduced by half, the situation is simplified greatly. It is 
observed that including the pion mass splitting modifies the binding energy by at most 0.5 MeV. We suggest that the 
2 ++ , ++ and h B*B* molecules should exist, they would be narrow states, and they dominantly decay into 
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APPENDIX A: THE ONE BOSON EXCHANGE POTENTIAL FOR THE D*D*{B*B*) STATES 

INCLUDING THE ISOSPIN MASS SPLITTING 
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Va~(r) = [Vc(r)+V s (r)] 
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In the above expressions, the parameters m n , and /i p could take two different sets of values. For ]J*°]J* — > jj*Ojj*o 
(B*°B*° -> B*°B*°) and D*+D*~ -► D*+D*~ (B*+B*- -> B*+B*-), we should choose TO ff = m^o, ^ = m^o and 
Whereas for the processes D*°D* -> D*+ D*- (B*°B*° -> B*+B*-) and D*+D*- -> (B*+B*- -> 

-(m £) »+-TO £ ,»o) 2 ] 1 / 2 = m,r±) and /Up = [m 2 -(m I) »+-TO £ ,,o) 2 ] :L/2 



B*°B*°), we should take = nv±,/i OT 
(/ip = m p ) respectively. 



APPENDIX B: NUMERICAL RESULTS FOR THE D*D* STATES 



16 



A(MeV) 


M(MeV) 


r rms (fm) 


P!f :PS':Pr :P£-(%) 


930 


4013.80 


8.24 


95.69:1.80:0.50:2.01 


940 


4011.72 


2.39 


70.69:11.01:7.43:10.88 


950 


4004.40 


1.36 


39.99:22.32:16.38:21.31 


960 


3990.26 


1.03 


25.21:29.74:16.55:28.50 


970 


3968.64 


0.86 


17.70:34.73:14.01:33.57 


980 


3938.98 


0.76 


13.19:38.20:11.41:37.19 


990 


3900.83 


0.68 


10.22:40.68:9.27:39.83 


all couplings 


are reduced by half except q^nn 


A(MeV) 


M(MeV) 


r rms (fm) 


P|?:PS':Pj-:PS-(%) 


1100 


4011.77 


2.36 


64.96:12.17:10.77:12.10 


1130 


3998.69 


1.15 


30.86:26.00:18.08:25.06 


1160 


3971.28 


0.85 


18.83:33.55:15.04:32.59 


1190 


3927.61 


0.70 


13.03:38.09:11.59:37.29 


1220 


3866.06 


0.61 


9.66:40.98:9.01:40.35 



TABLE II: The predictions for the static properties of the J — ++ D*D* hadronic molecule, where M denotes the mass, 
rms is the root of mean square radius, Ps and Pd represent the S state and D state probabilities respectively. 



A(MeV) 


M(MeV) 


r rms (fm) 


P P U : P£-(%) 


950 


4011.33 


1.76 


58.98:41.02 


970 


3995.72 


1.13 


53.55:46.45 


990 


3969.26 


0.89 


51.99:48.01 


1010 


3930.16 


0.74 


51.28:48.72 


1030 


3876.73 


0.64 


50.88:49.12 


all couplings except g 


irjviv are 


reduced by half 


A(MeV) 


M(MeV) 


r rm s (fm) 


P P U : Pp~(%) 


1120 


4012.74 


2.05 


60.78:39.22 


1150 


4004.34 


1.30 


54.77:45.23 


1200 


3980.14 


0.94 


52.20:47.80 


1250 


3940.92 


0.75 


51.28:48.72 


1300 


3883.95 


0.63 


50.83:49.17 



TABLE III: The predictions about the mass, the root of mean square radius(rms) and the probabilities of the different compo- 
nents for the 0~ + D*D* molecule. 



A(MeV) 


M(MeV) 


r rms (fm) 


pUU . p UU . pUU . p+- . p + - . p+- 


(%) 


920 


4006.82 


1.35 


4.68:50.36:0.42:3.59:40.55:0.40 




930 


3997.17 


1.11 


3.79:49.19:0.47:3.21:42.89:0.45 




940 


3984.57 


0.97 


3.21:48.71:0.51:2.86:44.19:0.50 




950 


3968.85 


0.87 


2.80:48.47:0.57:2.57:45.04:0.55 




960 


3949.86 


0.79 


2.49:48.33:0.62:2.33:45.61:0.61 




970 


3927.45 


0.73 


2.25:48.21:0.70:2.13:46.02:0.69 




980 


3901.44 


0.68 


2.06:48.11:0.79:1.97:46.29:0.78 






all couplings except g 


ttNN are reduced by half 




A(MeV) 


M(MeV) 


r rms (frn) 


p00 . pOO . pOU . p +- . p+- . p+- 
r P() • *P2 ■ r F • r P() ■ r P2 • r F 


(%) 


1050 


4007.15 


1.34 


5.36:49.41:0.49:4.16:40.10:0.48 




1080 


3991.87 


1.01 


4.10:48.03:0.59:3.60:43.11:0.58 




1110 


3969.82 


0.84 


3.38:47.66:0.67:3.12:44.51:0.66 




1140 


3940.40 


0.73 


2.90:47.52:0.77:2.75:45.31:0.75 




1170 


3903.00 


0.65 


2.56:47.43:0.88:2.46:45.79:0.87 





TABLE IV: The predictions about the mass, the root of mean square radius(rms) and the probabilities of the different compo- 
nents for the 1 _ ~ D*D* molecule, where P P q and Pp° denote the x Pi and 5 Pi D*°E)* states probabilities respectively. 
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A(MeV) 


M(MeV) 


r rms (fm 


\ pUU . pUU . pUU . pUU . pH — . p+- . pH — . p + - 


(%) 


860 


4010.82 


1.96 


0.39:67.73:1.95:0.01:0.38:27.59:1.96:0.01 




890 


4007.14 


1.49 


0.51:60.38:2.68:0.01:0.49:33.27:2.64:0.01 




920 


4001.88 


1.26 


0.64:55.35:3.60:0.02:0.61:36.24:3.51:0.02 




950 


3994.34 


1.10 


0.79:51.20:4.86:0.04:0.76:37.58:4.72:0.04 




dsn 
you 




u.yo 


n QS-/17 iQ-fi fi^-n no-n fi^-fi /i^-n no. 
u.yo.^i < . iy .o.uo.u.uy .u.y^.o < .oo.o.^o.u.uy 




1010 


3966.58 


0.88 


1.21:42.85:9.11:0.22:1.17:36.37:8.85:0.22 




1040 


3940.41 


0.79 


1.51:37.53:12.49:0.69:1.47:33.44:12.18:0.68 




1070 


3894 04 


0.69 


2.01:29.16:16.55:3.55:1.96:27.02:16.22:3.53 






all 


coupling 


s except Qtvnn are reduced by half 




A(MeV) 


M(MeV) 


r rms (fm 


\ puu . pUU . pUU . pUU . pH . p+- . pH . p+- 

I r D0 • r S ■ r D2 • r G ■ r D0 ■ r S 1 r D2 • r G 


(%) 


970 


4009.21 


1.66 


0.54:61.73:2.69:0.01:0.53:31.79:2.69:0.01 




i nnn 
1000 


A{\(\Ci OA 

4U06.34 


1 A 1 

1.41 


n co.c7 rn.'i icn no.n fin .o /i on .o itc.n no 
0.03:5 ( .00:3.18:0. 02:0. 60:34. yU:3. 15:0.02 




1100 


3991.66 


1.01 


0.89:49.31:5.04:0.05:0.87:38.86:4.94:0.05 




1200 


3965.34 


0.81 


1.21:43.80:7.63:0.14:1.18:38.43:7.48:0.14 




1300 


3918.43 


0.68 


1.59:38.25:11.19:0.47:1.56:35.48:11.01:0.46 




1350 


3882.03 


0.62 


1.82:34.98:13.32:0.95:1.79:33.04:13.15:0.95 





TABLE V: The predictions about the mass, the root of mean square radius(rms) and the probabilities of the different components 
for the 2 ++ D*D* molecule. 



